Mediated Entanglement And Correlations In A Star Network Of Interacting Spins 
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We investigate analytically a star network of spins, in which all spins interact exclusively with 
a central spin through Heisenberg XX couplings of equal strength. We find that the central spin 
correlates and entangles the other spins at zero temperature to a degree that depends on the total 
number of spins. Surprisingly, the entanglement depends on the evenness or oddness of this number 
and some correlations are substantial even for an infinite collection of spins. We show how symmetric 
multi-party states for optimal sharing and splitting of entanglement can be obtained in this system 
using a magnetic field. 
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For a long time spin correlations in ID chains and 
higher dimensional lattices of interacting spins have been 
a subject of extensive interest Recently, the same 

systems have been studied from the point of view of truly 
quantum correlations or entanglement [^-p^. However, 
lattices of various dimensions are not the only physical 
systems whose fabrication is possible with current tech- 
nology. In particular, with the advent of quantum com- 
putation, various technologies have evolved which can 
make any member of an array of qubits (systems isomor- 
phic to spin-1/2) controUably interact with any other 
member |l3|-[l5|. It thus becomes possible to visualize 
structures of interacting spins which do not fall into the 
category of lattices in various dimensions. One very sim- 
ple structure that one can imagine, is a spin-star, as op- 
posed to the extensively studied spin- chains. In such a 
spin star, there is a preferred spin, which we call the 
central spin which interacts with all the other spins. All 
the non-central spins (which we will call the outer spins), 
on the other hand, do not directly interact amoiig them- 
selves. The structure is clearly depicted in Fig.|[ de- 
picts the central spin. The spins 1 — 5 interact only with 
the central spin and not with each other. The architec- 
ture is analogous to the star distribution networks used 
in communications. To our knowledge, not just entangle- 
ment and correlations, but also the statistical mechanics 
of such a structure remains unexplored. 

The star configuration with couplings of equal strength 
has many symmetry properties due to its invariance un- 
der the exchange of any two outer spins and we solve it 
exactly in the case of an Heisenberg XX interaction. The 
XX model was intensively investigated for spin chains by 
Lieb, Schultz and Mattis has been realized in recent 
years as an effective Hamiltonian in some systems . 
We find that the central spin mediates correlations and 
entanglement between the outer spins at zero tempera- 
ture to a degree that depends on the total number of 
spins in the spin-star. As expected, the entanglement 
goes down on average with the increase of the total num- 
ber of spins. But it shows curious oscillations with the 
evenness or oddness of this number. This is surprising, 
because in a star network, addition of an extra outer spin 
in the network is only expected to make it harder for the 
central spin to mediate entanglement. For the same rea- 
son, it is also surprising that we find that irrespective of 



the number of the outer spins, the solitary central spin 
is capable of imposing substantial spin ordering (corre- 
lation function > 1/2) in the X and Y directions. We 
also show that we can apply a magnetic field to our sys- 
tem to obtain multi-party states for optimal symmetric 
splitting and optimal symmetric sharing |l^ of en- 
tanglement as the ground state and as a simple derivative 
of the ground state respectively. 
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FIG. 1. This figure depicts the star configuration of spins. 
The spin labelled is the central spin, which interacts by the 
XX Heisenberg interaction with spins 1-5 around it. 
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The Hamiltonian which describes our system is given 
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where the summation over 'outer' refers to the outer 



spms, fJi; 



and Uiy denote the Ux 



and Gy Pauli operators 
for the ith outer spin and uqx and a^y denote the ax and 
Uy Pauli operators for the central spin. It can be shown 
that Jx = (1/2) I^outer'^i^;,^!/ = (1/2) Pouter 
Jz — (1/2) X^ring '^'^^ obey the standard angular momen- 
tum commutation relations (we have taken h, = \). This 
implies that the outer spins collectively behave as a single 
spin with spin operator J = iJx + jJy + kJz- It can be 
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shown that = + + commutes with H. It will 

J, y z 

also help to define the total angular momentum opera- 
tor F = (l/2)a-o + J, where cro = Ictox + jo'oy + kaoz 
and it can be shown that the z-componcnt, obeys 
[H^Fz] — 0, [J'^^Fz] — 0. Therefore simultaneous eigen- 
states of H , and F^ can be constructed. 

It is convenient to recast the Hamiltonian in Eq.|l| us- 
ing the raising and lowering operators a± = {jJx ± iay) 
and J± = (l/2)Eouter^± as 



H = J{<To+J- +cra-J+) 



(2) 



The above Hamiltonian thus represents a resonant in- 
teraction between a spin- (1/2) and a higher spin (with 
operator J) system. Inspired by the above form of H 
(in particular, its similarities with the Jaynes-Cummings 
model of quantum optics we conjecture the follow- 

ing state as the general form of its eigenstates 
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(|0) b',m)±|l) |j,m-l)) 



(3) 



where the first ket in each term denotes the central spin 
(|0) and |1) stand for the | — 1/2) and |l/2) spin states 
of the central spin), and the second ket is an eigenstate 
of J^. j is the quantum number associated with eigen- 
states of (eigenvalue of is j{j + 1)), and m is the 
quantum number for Jz . The way \j, m) is paired with 
central spin state |0) and |j, m— 1) with |1) in Eq.(||) 
means that the state is an eigenstate of Fz with eigen- 
value TO — 1/2. Eq.|| is valid for m = j to m — —j + 1. 
There are also two additional states where only one of 
the terms exist: |1) \ j,j), because |0) \ j,j + 1) does not 
exist, and similarly |0) 

We now seek to prove that states of the form in Eq. (|^) 
are indeed eigenstates by applying the Hamiltonian in 
Eq.^ to them. Applying H gives 

^(1/ V2) (1 1) J_ \j, to) ± |0) J+ I j, TO - 1)) 

= Jii/V2) (|i) + + \j, 

±|0) ^{j - (to - l))(j + (to - 1) + 1) |i,TO)) 

= ±Jy/{j + ni,){j -m + 1) 

x(1/\/2)(|0) b-,TO)±|l)|j,TO-l)), 



where the standard relations for J± |j, to) have been used 
in the second step. This confirms that states in Eq.|| are 
eigenstates, with energy eigenvalues 



E ^ ±Jy/{j + m){j - m + I) 



(4) 



Next we prove that these states account for all possible 
eigenstates, so that we have found a complete eigen-basis 
for H. 

The total number of eigenstates of the form Eq.H is 4j 
for a given value of j (as to runs from j to —j + 1 and 



there are two states for each to due to the ± in Eq.^) . In 
addition, there are the two states where one of the terms 
in Eqjl was missing. In total therefore, for a given value 
of j, there are 4j -I- 2 possible eigenstates. If we can now 
enumerate the possible values of j (which is the label for 
the total angular momentum of the outer spins), we can 
find out the number of eigenstates we have been able to 
account for. 

Suppose there are N outer spins. Let us first consider 
the different ways of generating different to values. Then 
there are ^Cq = 1 ways of having all these spins aligned 
in the same direction, giving to = N/2. If all the outer 
spins but one are aligned, then to = (A^ — l)/2 — 1/2, 
and there will be ^Ci ways of getting this value of to. 
For the general case of N — r spins aligned and r spins 
anti- aligned to them, then to = {N — r)/2 — r/2 and there 
are C^. ways of getting this value of to. 

We now have to group the above to values under suit- 
able j values. Consider the case where to = N/2. As 
there is only one arrangement of spins giving rise to this, 
there should be only one value of j = A^/2 to account for 
this. In other words, if there was more than one allowed 
j ~ N/2 value for the collection of spins, there would 
be more than one to = A^/2, which we know, is not the 
case. Next consider to = (A^ — 2).i. There are ^Ci such 
ways of obtaining this value of to. One of these ways 
will arise from when j = N.^, which means there must 
be ^Ci Co of j = (AT - 2).i to account for the re- 
maining ways of getting this value of to. This procedure 
for determining possible j continues until all to are ac- 
counted for. In general, there cLr6 ^ Cf — ^ C*r — 1 

ways of 

obtaining j — (A — 2r).i, with allowed values of r rang- 
ing from r = to r = A^/2 if A^ is even, or r = (A^ — l)/2 
if A is odd. 

Given these values and j, and that each value of j can 
produce 4j -I- 2 eigenstates, then a summation expressing 
the total number of states we have accounted for is 



^ (^a Cr- 



r=0 



(Af-2r).i 
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where x = N/2 if A^ is even, or {N- 1) /2 if N is odd. The 
result of the summation is determined by observing that 
successive terms in summation cancel out some of the 
preceding terms. The result is independent of whether A^ 
is even or odd, and is 2^+^. This means the total number 
of eigenstates of the type we identified in Eq.|| is 2^+-'^, 
which is equal to the dimension of the Hilbert space of 
all the A^ outer spins plus the central spin. Therefore the 
states in Eq.^ (when combined with states of the type 
|1) |j, j) and |0) |j, —j)) describe all the eigenstates of H. 

The ground state is the state with the lowest energy. 
To identify it we will first have to assume a sign of J7, 
which we take to be positive. All our results about en- 
tanglement and correlations will be exactly the same for 
negative J. In Eq.^ the lowest energy is obtained when j 
has its maximum possible value and to has its minimum 
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possible value. For the case A'' odd, the lowest energy is 
when m = i i.e. the eigenstate 



(l/V2)(|0)|7V/2,l/2>-|l)|Ar/2,-l/2)) 



(6) 



and if N is even, then in fact the ground state is degener- 
ate because there are two states with the lowest possible 
energy, when m = or m = 1. These are 
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(|0> \N/2,0)-\l) \N/2,-l)) 



V2 



|7V/2,l)-|l)|iV/2,0)) 



(7) 



The reason for the above difference between N even and 
N odd is that the two cases lead to an integral and half 
integral value of j respectively. When j is half integral, 
m = ±1/2 is allowed and gives an unique ground state. 
For j integral, the 0, —1 and 1, form two distinct j, m 
pairs to combine with the central spin- 1/2 particle to give 
two degenerate ground states. 

To compute entanglement and correlations, it is useful 
to have expressions in terms of the states of the individ- 
ual outer spins for these ground states. Let |0) and |1) 
stand for the | — 1/2) and |l/2) spin states of any outer 
spin. For N odd, the state \N/2, 1/2) is an equal super- 
position of all states with {N -I- 1)/2 ones and {N — l)/2 
zeros with no relative phase between them. The state 
|7V/2, -1/2) is the same type of state with {N - l)/2 
ones and {N + l)/2 zeros. For example, for iV = 3, 



13/2,1/2) 
13/2,-1/2) 



1 

V3 
1 
V3 



[|011) 
[|100) 



1101) 
1010) 



1110)] 
1001)] 
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There are similar expressions for the ground state for N 
even. The |A^/2,0) state is an equal superposition of all 
states with an equal number of zeros and ones, with no 
relative phase between the superposed states. \N/2, ±1) 
is the same type of state with A^/2±l ones and rest zeros. 

Given the ground state, we are able to calculate the 
entanglement between any two outer spins in this state 
{i.e. at zero temperature). The symmetry of the problem 
implies that the entanglement will be the same between 
any two outer spins. Since H and Fz commute, the form 
of the reduced density matrix for any two spins is B 



/V 

w z 

z X 

VO 2/ 



in the standard basis {|00) , |01) , |10) , |11)}. For such 
density matrices, a measure of entanglement called the 
concurrence is given by Q] 

C = 2max{|z| - 0} 

For N odd, the concurrence comes out as 



C = 2max{^,0} = l 



For the case of iV even, where there are two ground states, 
a similar procedure is followed, except that the reduced 
density matrix is now described as an equal mixture of 
the two states. This gives the concurrence as 



C — 2 maxj — 



-2N 2N{N-l) 



1 



N N{N-l) 



Thus the entanglement goes to zero as iV — > oo, which is 
expected, as the entanglement is mediated by the central 
spin. The total entangling capability (and thereby me- 
diating capability) of the central spin is divided among 
a larger and larger number of outer spins as N becomes 
larger. However, on going from an even N to an odd 

-|- 1 number of outer spins, the concurrence rises from 
1/N - 1/{N^ - N) to l/(iV + 1). Thus on increasing 
N there are curious oscillations in concurrence of am- 
plitude 2/{N{N - l)iN + 1)}. These oscillations, due 
to the different expressions for concurrence in the cases 
of even and odd N, is due to mixedness of the zero 
temperature density matrix for even N. On applica- 
tion of a magnetic field in the +Z direction, the state 
(|0) \N/2,0) - |1) |iV/2,-l)) becomes the non-degenerate 
ground state for even TV and the oscillations in entangle- 
ment disappear. 

There are also interesting results to be noted for cor- 
relation functions in the star network. The {crizcr2z) cor- 
relations follow the same pattern as the entanglement, 
but 



i + ^foroddiV 

1 1 1 

2 2iV ~ 2A^(A^- 1) 



for even N . 



What is particularly surprising above, is the non- 
vanishing nature of the correlations in the large N limit. 
The solitary central spin is capable to imposing spin or- 
der in the X direction (so that {(Jix(J2x) = 1/2), even 
when there are an infinite number of outer spins to order. 
The same result holds for {criya2y). Thus our system pro- 
vides an effective way of imposing order simultaneously in 
the X and Y directions for an infinite collection of spins. 
This result also highlights a crucial difference between en- 
tanglement and correlations: while a finite dimensional 
quantum system cannot be individually entangled to each 
member of an infinite collection of systems, it can indeed 
be correlated individually to each of them. 

We now show that the application of a magnetic field 
allows us to change the ground state to 



l«) = ^(|0)|7V/2,-iV/2+l)-|l)|7V/2,- 
= ^(|0){|000...1)}-|1)|000...0)) 



-N/2)) 



(9) 
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where {|000 • ■ ■ 1)} is a normalized state that is an equal 
superposition of all states with only one |1) with no rel- 
ative phase between the superposed components. This 
state has the special significance that the concurrence 
between the central spin and each of the outer spins is 
1 / ViV, which is the maximum consistent with symmet- 
ric splitting of the total entanglement of one qubit with 
a collection of N qubits among the N qubits . To our 
knowledge, this is the first identification of the canonical 
state I a) as the ground state of an interacting spin sys- 
tem. Once the ground state |a) is generated, if the cen- 
tral spin is measured and found to be in the state |0), the 
rest of the spins are projected onto the state {|000 ... 1)}. 
The central spin can now be removed to make the state 
dynamically steady (except for decoherence and sponta- 
neous decay effects). This state has the property that 
the concurrence between any two spins is 2/N, which is 
the maximum possible entanglement in a collection of N 
spins in which all pairs of spins are equally entangled . 
While in Ref. ||^, it was only conjectured that a state of 
the type {|000 ■..!)} could be made the ground state of 
the isotropic closed Heisenberg chain using a magnetic 
field, here we will rigorously prove the preparation of \a) 
and thereby {|000...1)}. 

To show that it is possible to make |a) the ground 
state, we consider the energy eigenvalues in a uniform 
magnetic field B in the +Z direction (in which the eigen- 
states remain unchanged) 



1, 



rn + 1) + (m — 



(10) 



When B becomes so high that the second term in E 
dominates, the relative ordering of the energy levels will 
be determined purely by m, with the ground state be- 
ing state with m = i.e. \(3) = |0) |000 . . . 0). It 
is straightforward to show that has an energy lower 
than I a) for B > JypN . Before |/3) becomes the ground 
state, there is a range of B in which |a) is the ground 
state. To prove this, wc will have to show that the energy 
of I a) will be less than that of all other states in a certain 
range. The value of B for which \a) has the same energy 
as a general state described by j and m is given by 



B = 



J (^[N - v/(j+m)(j-m+l)) 



(11) 



maxi- 
This 
2 for 



B attains its largest value, when j is a 
mum and m has it's most negative value, 
happens when j ~ and m = — -y 4 

which B = Jy^(^y/2{1 - 1/iV) - l) < J^/N. 

Therefore, for a magnetic field with a range 

J^/N ( J2Tr^^l7N) - l] < B < JVN, the state la) 



N. 
of 

/7V(^V2(l-l/iV) ' ~ 

is the ground state. 

In this letter we have introduced and solved a spin- 
star, an architecture of interacting spins which cannot be 



classified as a lattice in any dimension. It is physically re- 
alizable in various arrays of qubits designed for quantum 
computation. Testing for the entanglement and corre- 
lations predicted here would serve as a benchmark test 
for the functioning of arrays of qubits. If spin-spin inter- 
actions can be extended to long distances, the spin-star 
could be used for distribution of entanglement through 
the multiparty states for optimal symmetric sharing and 
splitting of entanglement. Exploration of the full statis- 
tical mechanics of a spin-star, would be interesting future 
work. 

AH thanks the UK EPSRC (Engineering and Physical 
Sciences Research Council) for financial support. 



[2: 

[4: 
[6: 

[7 
[8 

[9: 
[10: 



[11 

[12: 
[is: 



[14 

[is: 
[16: 



[17: 
[is: 

[19 



H. A. Bethe, Z. Physik 71, 205 (1931); J. Bonner and M. 
E. Fisher, Phys. Rev. 135, A640 (1964). 
B. Kaufman and L. Onsager, Phys. Rev. 76, 1244 (1949). 
E. H. Lieb, T.D. Schultz, D.C. Mattis, Annals of Phys. 
(N.Y.) 16, 407 (1961). 

K. M. O'Conner and W. K. Wootters, Phys. Rev. A 63, 
052302 (2001). 

M. A. Nielsen, Ph .D thesis, University of Mexico, 1998, 
quant-ph/001103e . 

M. C. Arnesen, S. Bose, and V. Vedral, Phys. Rev. Lett. 
87, 017901 (2001). 

X. Wang, Phys. Rev. A 64, 012313 (2001); X. Wang, 
Phys. Lett. A 281 101 (2001); X. Wang, H. Fu and A. I. 
Solomon, J. Phys A-Math Gen. 35, 4293 (2002). 
D. Gunlycke, V. Kendon, V. Vedral and S. Bose, Phys. 
Rev. A 64, 042302 (2001). 

G. L. Kamta and A. F. Starace, Phys. Rev. Lett. 88, 
107901 (2002). 

A. Osterloh, L. Amico, G. Falci and R. Fazio, Nature 



416, 608 (20 02); T. J. Osborne and M. A. Nielsen, quant 



ph/0202162 



X. Wang, quantph /0203141; X. Wang and P. Zanardi, 
quant-ph/020210^. 



I. Bose and E. Chattopadhyay, 3ond-mat/0208011 



J. I. Cirac and P. Zoller, Nature 404, 579 (2000); B. E. 
Kane, Nature 393, 133 (1998); Y. Makhlin and G. Schon, 
Nature 398, 305 (1999). 

A. Imamoglu et. al, Phys. Rev. Lett. 83, 4204 (1999). 
S. B. Zheng and G. C. Guo, Phys. Rev. Lett. 85, 2392 

(2000) . 

V. Coffman, J. Kundu and W. K. Wootters, Phys. Rev. 
A 61, 052306 (2000); W. Diir, Phys. Rev. A 63, 020303 

(2001) . 

M. Koashi, V. Buzek, and N. Imoto, Phys. Rev. A 62, 
050302 (2000). 

B. W. Shore and P. L. Knight, J. Mod. Opt. 40, 1195 
(1993). 

W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 



4 



